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Abstract 

Properties of mesons are investigated within the Nambu-Jona-Lasinio model. 
We include meson-loop corrections, which are generated via a systematic 1/Nc- 
expansion in next-to-leading order. We show that our scheme is consistent with 
chiral symmetry, in particular with the Goldstone theorem and the Gell-Mann Oakes 
Renner relation. The numerical part focuses on the pion and the p- meson sector. For 
the latter the l/A'^c-corrections are crucial in order to include the dominant p vrvr- 
decay channel, while the leading-order approximation only contains unphysical qq- 
decay channels. We show that a satisfactory description of the pion electromagnetic 
form factor can be obtained. Similarities and differences to hadronic models are 
discussed. 
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1 Introduction 



The understanding of hadron properties in the vacuum as well as in hot or dense matter 
is one of the central tasks of present-day nuclear physics. In principle, all properties 
of strongly interacting particles should be derived from QCD. However, at least in the 
low-energy regime, where perturbation theory is not applicable, this is presently limited 
to a rather small number of observables which can be studied on the lattice, while more 
complex processes have to be described within model calculations. 

So far the best descriptions of hadronic spectra, decays and scattering processes are 
obtained within phenomenological hadronic models. For instance the pion electromagnetic 
form factor in the time-like region can be reproduced rather well within a simple vector 
dominance model with a dressed p-meson which is constructed by coupling a bare p- 
meson to a two-pion intermediate state 0] . Models of this type have been successfully 
extended to investigate medium modifications of vector mesons and to calculate dilepton 
production rates in hot and dense hadronic matter P]. 

In this situation one might ask how the phenomenologially successful hadronic models 
emerge from the underlying quark structure and the symmetry properties of QCD. Since 
this question cannot be answered at present from first principles it has to be addressed 
within quark models. For light hadrons chiral symmetry and its spontaneous breaking 
in the physical vacuum through instantons plays the decisive role in describing the two- 
point correlators |Q with confinement being much less important. This feature is captured 
by the Nambu-Jona-Lasinio(NJL) model in which the four-fermion interactions can be 
viewed as being induced by instantons. 

The study of hadrons within the NJL model has of course a long history. In fact, 
mesons of various quantum numbers have already been discussed in the original papers 
by Nambu and Jona-Lasinio p and by many authors thereafter (for reviews see [|, |3, H). 
Most of these works correspond to a leading-order approximation in l/A^^c, the inverse 
number of colors. In this scheme quark masses are calculated in mean-field approximation 
and mesons are constructed as correlated quark- ant iquark states. With the appropriate 
choice of parameters chiral symmetry, which is an (approximate) symmetry of the model 
Lagrangian, is spontaneously broken in the vacuum and pions emerge as (nearly) massless 
Goldstone bosons. While this is clearly one of the successes of the model, the description 
of other mesons is more problematic. One reason is the fact that the NJL model does 
not confine quarks. As a consequence a meson can decay into free constituent quarks if 
its mass is larger than twice the constituent quark mass m. Hence, for a typical value 
of m ~ 300 MeV, the p-meson with a mass of 770 MeV, for instance, would be unstable 
against decay into quarks. On the other hand the physical decay channel of the p-meson 
into two pions is not included in the standard approximation. Hence, even if a large 
constituent quark mass is chosen in order to suppress the unphysical decays into quarks, 
one obtains a poor description of the p-meson propagator and related observables, like 
the pion electromagnetic form factor. 

This and other reasons have motivated several authors to go beyond the standard 
approximation scheme and to include mesonic fluctuations. In Ref. a quark- ant iquark 
p-meson is coupled via a quark triangle to a two-pion state. Also higher-order corrections 
to the quark self-energy and to the quark condensate have been investigated. 
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However, as the most important feature of the NJL model is chiral symmetry, one should 
use an approximation scheme which conserves the symmetry properties, to ensure the 
existence of massless Goldstone bosons. 

Two slightly different symmetry conserving approximation schemes are discussed in 
Refs. |T^ and [jl3|. The authors of Ref. use an expansion in 1/Nc in order to include 



mesonic fluctuations and calculate the changes of and of the quark condensate (V'V') 
a low-momentum expansion for the incoming state. In Ref. a correction term to the 
quark self-energy is included in the gap equation. In a perturbative scheme this term would 
be of order 1/Nc but, as the modified gap equation is solved self-consistently, arbitrary 
orders in l/N^ are generated. The authors find a consistent scheme to describe mesons 
and show the validity of the Goldstone theorem and the Goldberger-Treiman relation in 
that scheme. Based on this model various authors have investigated the effect of meson- 
loop corrections on the pion electromagnetic form factor |14] and on tt-tt scattering in the 
vector |]TB[ and the scalar channel |]TB[. However, since the numerical calculation of the 



multi-loop diagrams, which have to be evaluated, is rather involved, in these references 
the exact expressions are approximated by low-momentum expansions. 

In the present paper we calculate the p-meson self energy in a strict expansion up 
to next-to-leading order without any approximations. Within the same scheme we have 



recently studied the influence of mesonic fluctuations on the pion propagator [17|. This 



was mainly motivated by recent works by Kleinert and Van den Bossche ||T^, who claim 
that chiral symmetry is not spontaneously broken in the NJL model as a result of strong 
mesonic fluctuations. In Ref. [0 we argue that because of the non-renormalizability of 
the NJL model new divergences and hence new cutoff parameters emerge if one includes 



meson loops. Following Refs. [12| and |13] we regularize the meson loops by an inde- 



pendent cutoff parameter Am- The results are, of course, strongly dependent on this 
parameter. Whereas for moderate values of the pion properties change only quantita- 
tively strong instabilities are encountered for larger values of Aa/, which might be a hint 
for an instability of the spontaneously broken vacuum state. 

In Ref. [jl^ we restricted ourselves to calculate the pion mass m^, the pion decay 
constant and the quark condensate {ipip). Since these quantities can already be fitted 
without meson- loop corrections (i.e. A^ = 0) one has to look at other observables in 
order to fix Km- Obviously the spectral function of the p-meson is particularly suited, 
as it cannot be described realistically without taking into account pion loops. In the 
present article all parameters of the NJL model will be fixed by fitting p-meson properties 
simultaneously with f^^ and (ipip). The important result is that such a fit can indeed be 
achieved with a constituent quark mass which is large enough such that the unphysical 
gg-threshold opens above the p-meson peak. Since the constituent quark mass is not an 
independent input parameter this was not clear a priori. 

The paper is organized as follows. In Sec. |^ we present the scheme for describing 
mesons in next-to-leading order in l/N^. The consistency of this scheme with the Gold- 
stone theorem and with the Gell-Mann Oakes Renner relation will be shown in Sec. |^. 
In Sec. ^ we perform a low-momentum expansion of the effective meson vertices and dis- 
cuss the relation to hadronic models. The numerical results will be presented in Sec. |^. 
For several values of Am we fix a subset of the model parameters by fitting quantities in 
the pion sector. Then we determine the remaining parameters, in particular Am, in the 
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p-meson sector. Finally, conclusions are drawn in Sec. ^. 



2 The NJL model in leading order and next-to-leading 
order in l/N^ 

We consider the following NJL-type Lagrangian: 

C = - mo)V^ + gs + {^il^T^f\ - g. {{^l^'r^f + {^l^'l^r^f] ■ (1) 

Here ^/^ is a quark field with N j = 2 flavors and N^. colors, while gg and gy are dimensionful 
coupling constants of the order l/Nc- In order to establish the counting scheme, the 
number of colors has been treated as variable, but all numerical calculations will be done 
with the physical value, = 3. 

In most publications the NJL model has been treated in leading order in l/Nc- In 
terms of many-body theory this corresponds to a (Bogoliubov) Hartree approximation for 
the quark propagator and to a random phase approximation (RPA) for describing mesonic 
excitations. Diagrammatically this is shown in Fig. |l|. The selfconsistent solution of the 
Dyson equation shown in the upper part of Fig. [l| leads to a momentum independent 
quark self energy and therefore only gives a correction to the constituent quark mass: 



Usually one refers to this equation as the gap equation. For sufficiently large couplings gg 
it allows for a finite constituent quark mass m even in the chiral limit, i.e. for rriQ = 0. In 
the mean-field approximation this solution minimizes the ground-state energy. Since gg 
is of order l/N^. the constituent quark mass m, and hence the quark propagator S{p) = 
— m)~^ are of the order unity. 

Mesons are described via a Bethe-Salpeter equation, as shown in the lower part of 
Fig. |I[ Since this is all standard we only list the results, which are needed later on. First 
we define the quark-antiquark polarization functions 



Figure 1: (Upper part) The Dyson equation for the quark propagator in the Hartree ap- 
proximation (solid lines). The dashed lines denote the hare quark propagator. (Lower 
part) The Bethe-Salpeter equation for the meson propagator in the RPA (double line). 
Again, the single solid lines indicate quark propagators in the Hartree approximation. 





(3) 
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with M = a,7r,p,ai and = 1, = i'y^r^, r^*= = 'yi^r^ and T^l" = Yl^^^- Here "Tr" 
denotes a trace in color, flavor and Dirac space. Hm is diagrammatically shown in Fig. |^. 
Iterating the scalar (pseudoscalar) part of the four-fermion interaction one obtains for the 
sigma meson (pion): 

DM) = , r W ^f(g) ^ DM)5ah = . ~\ ^ , 5ab. (4) 



ab- 



Here a and h are isospin indices and we have used the notation (g) = n^(g) 5, 

In the vector channel this can be done in a similar way. Using the transverse structure 
of the polarization loop in the vector channel, 

m^^''\q) = n,(g) T^^-5ab ; T'^" = {-Q^" + ^) , (5) 



q2 



one obtains for the p-meson 



D^r\q) = D,{q) T^" Sab = ^^^^^ ' 

Analogously, the Oi can be constructed from the transverse part of the axial polarization 
function Ilaj. As discussed e.g. in Ref. U^'^ also contains a longitudinal part which 
contributes to the pion. Although there is no conceptional problem to include this mixing 
we will neglect it in the present paper in order to keep the structure of the model as 
simple as possible. 

It follows from Eqs. (Q) - (H) that the functions -Dm('?) are of order l/Nc. Their explicit 
forms are given in App. |B[ For simplicity we will call them "propagators", although 
strictly speaking, they should be interpreted as the product of a renormalized meson 
propagator with a squared quark- meson coupling constant. The latter is given by the 
inverse residue of the function DM{q), while the pole position determines the meson 
mass: 

n-l/ M n -2(0) dllMiq). , . 



We have used the superscript (0) to indicate that rn'l^p and Qf^j^^ are leading-order quan- 



titles in 1/iVc. One easily verifles that they are of order unity and l/^/Wl, respectively. 
With the help of the gap equation, Eq. (^, one can show that the pion is massless in the 
chiral limit, demonstrating the consistency of the scheme with chiral symmetry 0. 

We now turn to the corrections in next-to-leading order in l/N^. The correction terms 
to the quark self-energy are shown in Fig. 0. In these diagrams the quark lines and the 




a, 



Figure 2: The 1/ N^- corrections to the quark self-energy. 
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RPA (a) (b) (c) (d) 

Figure 3: Contributions to the mesonic polarization function in leading (RPA) and next- 
to-leading order in l/Nc- 

double lines correspond to quark propagators in the Hartree approximation (order unity) 
and to meson propagators in the RPA (order 1/Nc), respectively. Both diagrams are 
therefore of order l/Nc- 

The l/A^^c-corrected mesonic polarization diagrams read 

nM(g) = nM(g) + Yl ^n^^g) . (8) 

k=a,b,c,d 

The four correction terms SH^^j - SH^^ together with the leading-order term 11^ are shown 
in Fig. 1^. Again the lines in this figure correspond to Hartree quarks and RPA mesons. 
Besides the RPA meson propagators the main building blocks are quark triangle and box 
diagrams, which are shown in Fig. ^. The triangle diagrams entering into 611^^ and 511^'' 
can be interpreted as effective three- meson vertices. For external mesons Mi, M2 and M3 
they are given by 

/d'^k f 
j^-^^TT[TMiiS{k)TM2iS{k - p)TM,iS{k + q)] 

+TT[rM,tS{k - q)^M,^S{k + p)^M,^S{k)]] , (9) 



with the operators Tm as defined below Eq. (^. We have summed over both possible 
orientations of the quark loop. 

The quark box diagrams are effective four-meson vertices and are needed for the eval- 
uation of 5n^/ and ^n^v/ • If one again sums over both orientations of the quark loop they 
are given by 

—i^Mi,M2,M3,M4, (Pl! P2, Pa) 

d'^k 



TT[TM^iS{k)TM2'i'S{k-p2)TMaiS{k-p2-p3)TAn'i'S{k-^Pi)] 



{2nr 

-^TT[TMiiS{k - pi)TMjS{k + p2 + p3)TM3iS{k P2)TM2iS{k)]] . (10) 



The polarization diagrams SH^^j - 611^^^ contain several loops. However, with the help 
of the definitions given above they can be written in a relatively compact form: 

^nif(g) = ^ f rM,Mi,M2(g,p)^Afi(p)rAf,Mi,M2(-g, -p)^M2(-p-g) , 

^ ' M1M2 
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,P2,m\ / 



\p' = -p-q,M3 /'pi,M, 



^ p-i,M-i 



\ P4,M4 



Figure 4: (Left) The quark triangle vertex —iTMi,M2,M3il^P)- (Right) The quark box vertex 



—I 



2 J {27vy 



Ml 



(27r) 



TmiMiAp^ -p) DmAp) 



(11) 



Ml 



Note that for and 5Iif) one has to include a symmetry factor of 1/2, because 

otherwise the sum over the two orientations of the quark propagators, which is contained 
in the definitions of the quark triangle vertex (Eq. (^)) and of the quark box vertex 
(Eq. (pUD) would lead to double counting. Similarly in ^n^-* we had to correct for the fact 
that the exchange of Mi and M2 leads to identical diagrams. 

For the evaluation of Eq. (|TTp we have to proceed in two steps. In the first step we 
calculate the intermediate RPA meson-propagators. Simultaneously we can calculate the 
quark triangles and box diagrams. One is then left with a meson loop which has to be 
evaluated in a second step. 

The various sums are, in principle, over all quantum numbers of the intermediate 
mesons. Of course, many combinations vanish, e.g. because of isospin conservation. In 
fact, the expression for 511^^1 {q) should contain a sum over the quantum numbers of both 
intermediate mesons. However, one can easily verify that the meson which connects the 
two quark triangles has to be a sigma meson in order to give a non-vanishing contribution. 

As we have mentioned earlier, the main focus of the present paper is on the p-meson. In 
this case the most important contribution comes from the two-pion intermediate state in 
diagram 5n^"\ Other contributions to this diagram, i.e vrai, per, pp and aiOi intermediate 
states, are much less important since the corresponding decay channels open far above the 
p-meson mass and - in the NJL model - also above the unphysical two-quark threshold. 
Hence, from a purely phenomenological point of view, it should be sufficient to restrict 
the sums in Eq. ([TID to intermediate pions. However, in order to stay consistent with 
chiral symmetry, we have to include intermediate sigma mesons as well. On the other 
hand, vector- and axial-vector mesons can be neglected without violating chiral symmetry. 
Since this leads to an appreciable simplification of the numerics we have restricted the 
intermediate degrees of freedom to scalar and pseudoscalar mesons in the present paper. 
Of course, in order to describe a p-meson, we have to take vector couplings at the external 
vertices of the diagrams shown in Fig. R|. 
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In analogy to Eqs. (0) and the corrected meson propagators are given by 

= ^ , , . (12) 

1 - 2gMilM{q) 

with qm = Qs for a and vr and qm = Qv for P and ai. Since our scheme corresponds to a 
1 /Ai'c-expansion of the polarization diagrams and hence of the inverse meson propagators, 
the Dm contains terms of arbitrary orders in 1/iVc. 

The corrected meson masses are again defined by the pole positions of the propagators: 

D-^{q)U^=ml, = 0. (13) 

As we will show in the next section our scheme is consistent with the Goldstone theorem, 
i.e. in the chiral limit it leads to massless pions. Note, however, that because of its 
implicit definition itim contains terms of arbitrary orders in l/Nc-, although we start from 
a strict expansion of the inverse meson propagator up to next-to-leading order. This will 
be important in the context of the Gell-Mann Oakes Renner relation. 



Finally we should comment on Ref. |T2|, where the authors introduce 1/A^c-corrections 
in a slightly different way. In that scheme a correction term to the quark self-energy 
is self-consistently included in the gap equation. This leads to a modified quark mass 
as compared to the 'Hartree mass'. Mesons are described by iterating the polarization 
diagrams IIm, ^'^^m ■> ^'^^m ^^"^ ^'^^mi but using the modified quarks for constructing RPA 
mesons as well as quark triangle- and box diagrams. Since the l/A'^c-correction terms 
are iterated in the gap equation, all diagrams contain terms of arbitrary orders in l/iV^. 
However, if one performs a strict l/A^^c-expansion of the mesonic polarization diagrams 
one exactly recovers our scheme of describing mesons [|l^ . 



Both schemes, the one introduced in Ref. [0 and the strict l/A^^c-expansion, are 
consistent with chiral symmetry. In particular, they lead to massless pions in the chiral 



limit. However, because of the modified gap equation, in the scheme of Ref. |]T2[ the RPA 
pions are not massless in the chiral limit, but tachyonic. Therefore this model is not 
very well suited for calculating the p-meson self-energy, which is mainly determined by 
intermediate RPA pions. For this reason we prefer the strict l/A^^c-expansion scheme, in 
which both RPA- and l/A'^c-corrected pions are massless in the chiral limit. Still there 
remains the problem that ^ Tn^\ if we go away from the chiral limit, but at least for 
our final parameter set the difference is only about 10%. This will be discussed in more 
detail in Sec. |5|. 



3 Consistency with chiral symmetry 

Before discussing the numerical details we wish to show that the scheme introduced in 
the previous section is consistent with chiral symmetry. We begin with the Goldstone 
theorem and then show the consistency of the scheme with the Gell-Mann Oakes Renner 
relation. This is not an entirely academic exercise. Since most of the integrals which 
have to be evaluated are divergent and hence must be regularized one has to ensure that 
the various symmetry relations are not destroyed by the regularization. To this end, it is 
important to know how these relations formally emerge. 
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The validity of the Goldstone theorem in our scheme has aheady been proven by 
Dmitrasinovic et ah [0 and we only summarize the main steps which are relevant for our 
later discussion. One has to show that, in the chiral limit, the inverse pion propagator 
vanishes at zero momentum, 

2^,11^(0) = 1 for mo = 0. (14) 

As before we use the notation 11"* = Sab^n- Restricting the calculation to the chiral limit 
and to zero momentum simplifies the expressions considerably and Eq. ( p!^ can be proven 
analytically. Since the Goldstone theorem is fulfilled in leading order, i.e. 2gsIlT^{0) = 1 
for mo = 0, we only need to show that the contributions of the correction terms add up 
to zero: 

J2 5ni^\0) = for mo = 0. (15) 

k=a,b,c,d 

Let us begin with diagram 5Il^\ As mentioned above, we neglect the p and ai subspace 
for intermediate mesons. Then one can easily see that the external pion can only couple 
to a ira intermediate state. Evaluating the trace in Eq. (|^) for zero external momentum 
one gets for the corresponding triangle diagram: 

K'^,AO,p) = -5ab ^N^Nj 2m I{p) , (16) 

with a and h being isospin indices and the elementary integral 

^^^^ (27r)4(F-m2 + ^5)((A;+p)2-m2 + i£) ' ^^^^ 
Inserting this into Eq. (|ri|) we find 

5Il^^^\Q) = i5,b j -^,{m,NfI{p))Hm' D^ip) D^{p) . (18) 

Now the essential step is to realize that the product of the RPA sigma- and pion propa- 
gators can be converted into a difference [Q, 

to finally obtain 

5nWa.(o) ^ 4^^^^ I _^2/(p){d.(p) - D^ip)] . (20) 
The next two diagrams can be evaluated straightforwardly. One finds: 

'^"(0) = -5ab I -0^,{Da{p) (/(P) + /(O) - ip' - 4m2) K{p)) 

+D^ip) (Slip) + 3/(0) - 3/ Kip)) } , 
Sui^)^\0) = -6abAN^Ns j ^J{p){-D^{p)-D^{p)]. (21) 



The elementary integral K{p) is of the same type as the integral I{p) and is defined in 
App. 0. Finally we have to calculate 6Il^\o). According to Eq. (|ll|), it can be written 
in the form 

^^""(0) = -zC^,,,(0,0)D.(0)A, (22) 

with 

^{ D^ip) (2 lip) + /(O) - (/ - 4m') K{p)) 

+D^ip) ( 3/(0) - 3p' K{p) ) } . (23) 

Evaluating D„{0) in the chiral limit and comparing the result with Eq. ([T^) one finds that 
the product of the first two factors in Eq. (p^) is simply 6ab/m, i.e. one gets 

5nWafe(0)=5,,-. (24) 

m 

With these results one can easily check that Eq. ([15|) indeed holds in our scheme. 

Of course, most of the integrals we have to deal with are divergent and have to be 
regularized. Therefore one has to make sure that all steps which lead to Eq. ([I5| ) remain 
valid in the regularized model. One important observation is that the cancellations occur 
already on the level of the p-integrand, i.e. before performing the meson-loop integral. 
This means that there is no restriction on the regularization of this loop. We also do 
not need to perform the various quark loop integrals explicitly but can make use of 
several relations between them. For instance, in order to arrive at Eq. ( p4|) we need the 
similar structure of quark triangle r7r,7r,o-(0, 0) and the inverse RPA propagator Du{0)~^. 
Therefore all quark loops, i.e. RPA polarizations, triangles and box diagrams should 
be consistently regularized within the same scheme, whereas the meson loops can be 
regularized independently. 

Going away from the chiral limit the pion gets a finite mass. To lowest order in the 
current quark mass it is given by the Gell-Mann Oakes Renner (GOR) relation, 

mlf^ = -mo{^^) . (25) 

Our scheme is consistent with this relation up to next-to- leading order in l/Nc- Expanding 
both sides of Eq. (^5|) we get 

Here rrh!'^^ denotes the leading-order contribution to the squared pion mass, while is 
the next-to-leading order contribution. For and (ipip) we introduce similar notations. 
Since the GOR relation holds only in lowest order in mo, Eq. (^6]) corresponds to a double 
expansion: has to be calculated in linear order in mo, and {ijjip) in the chiral limit. 
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The quark condensate is given by 



m = -^j ^J^Sip) , (27) 

with the l/Ai'c-corrected quark propagator 

S{p) = S{p) + + (28) 

Here S{p) is the Hartree propag ator and S^"^) and S(^) are the self energy corrections 
shown in Fig. ^. Since we are interested in a strict 1/Nc expansion, it should not be 
iterated. The r.h.s. of Eq. ( p7|) can be easily evaluated and one obtains 



(«.)'«)+ MM = - ^f^. (29) 

with A as defined in Eq. (p^). 

The pion decay constant is calculated from the one-pion to vacuum axial vector ma- 
trix element. Basically this corresponds to evaluating the mesonic polarization diagrams, 
Fig. 1^, coupled to an external axial current and to a pion. This leads to expressions 
similar to Eqs. (0) and (0), but with one external vertex equal to 7'^75^, corresponding 
to the axial current, and the second external vertex equal to gnqqilbT^, corresponding to 
the pion. Here the l/A^^c-corrected pion-quark coupling constant is defined as 



„-2 _ -2(0) , r -2 _ dU^iQ) I /oQ^ 
ynqq " Unqq + '^dnqq " \<l^=mi ) {^^ ) 



analogously to Eq. (^. Now we take the divergence of the axial current and then use the 
axial Ward identity 

75 1^ = 2m75 +j5S-\k + p) +S-\k)j, (31) 



to simplify the expressions [|T2|. One finds: 

/ fl^(g)-n^(0) H.(g) - H.(0) ^ ,^..\ 
/. = 9.qq [ rn + DM A ) • (32) 

In the chiral limit, = — > 0, Eqs. (|^) and (^) can be employed to replace the 
difference ratios on the r.h.s. by pion-quark coupling constants. When we square this 
result and only keep the leading order and the next-to-leading order in 1/Nc to finally 
obtain: 

fT^ + 5fl = m^g^qT + {^'^9-^, + 2mZ).(0)A^;,f)) . (33) 

For the pion mass we start from Eqs. (|12D and (|13]) and expand the inverse pion 
propagator around = 0: 

1 - 2g,UM - 2(7,(^n^(g))| ml + 0{ml) = 0. (34) 
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To find in lowest non- vanishing order in mo we have to expand 1 — 2gsIlTT{0) up to 
hnear order in mg, while the derivative has to be calculated in the chiral limit, where it 
can be identified with the inverse squared pion-quark coupling constant, Eq. (pO]). The 
result can be written in the form 

.4 = ^f^(l-^)+0«. (35) 

Finally one has to expand this equation in powers of l/Nc- This amounts to expanding 
g'^qq, which is the only term in Eq. (^) which is not of a definite order in l/iV^. One gets: 



^2(0) , . 2 _ _ mn — ^("o^WAo-S + Djfi)A\ .of^X 

+ dm^ - mo ruo ^g^^^ dg^^^ + ^ j. [6b) 

It can be seen immediately that the leading-order term is exactly equal to —mQ^ipip)'^^^ / f'i'^^\ 
as required by the GOR relation. Moreover, combining Eqs. ( p9D for mo = 0, (^) and (|36| ) 
one finds that the GOR relation in next-to-leading order, Eq. (p6D, holds in our scheme. 

However, one should emphasize that this result is obtained by a strict 1/iVc-expansion 
of the various properties which enter into the GOR relation and of the GOR relation 
itself. If one takes and m^ as they result from Eqs. (|3^) and ( pSf ) and inserts them into 



the l.h.s. of Eq. ( pS] ) one will in general find deviations from the r.h.s. which are due to 
higher-order terms in l/N^. In this sense one can take the violation of the GOR relation 



as a measure for the importance of these higher-order terms [1/ 



4 Relation to hadronic models 

This section discusses an approximation to our scheme which points out the relation 
to hadronic models. In order to suppress the quark effects in the present model, it 
is suggestive to assume that the constituent quark mass is very large as compared to 
the relevant meson momenta. This assumption leads to an effective zero momentum 
approximation, i.e. the quark vertices are taken at zero incoming momentum ("static 
limit"). In order to preserve chiral symmetry we then have to approximate the RPA- 
meson propagators consistently. This amounts to replacing the function J(p) in the RPA 
polarization functions by /(O). The latter can be related to the leading-order pion decay 
constant in this approximation. The RPA-meson propagators then take the form of free 
boson propagators: 

2(0) 

D'Sd) = . (37) 

— my 

with 



2gj:'"> AgJ.. 

and 



(38) 



2(0) ^ 2(0) ^ 2 2(0) ^ 2 2(0) ^ rn^ 

tl-n-qq iJaqq o^pqq o^aiqq ^2(0) ' ^ ' 
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Next we expand the quark triangles and box diagrams to first order in the external 
momenta. In this way one obtains a hadronic model with effective meson-meson coupling 
constants. Let us, for example, look at the p-meson self-energy which is generated as an 
approximation to the l/A^^c-corrected NJL model. As the ppa- and the ppaa-vertex vanish 
in this approximation, we are left with a pion loop diagram, which is generated from the 
diagram shown in Fig. |^(a), and a pion tadpole diagram, which arises from the sum of 
the diagrams in Fig. |^(b) and (c). These are exactly the diagrams which are calculated 



in standard hadronic descriptions for the p-meson in vacuum, e.g. 0, pO |. 

In hadronic models one is also interested in the coupling of a bare p-meson to a photon, 
which is e.g. needed to calculate the pion electromagnetic form factor via vector-meson 
dominance. In the NJL model the bare p-meson corresponds to the RPA meson and 
its coupling to a photon is basically given by the RPA polarization loop in the vector 
channel. If we now perform the same low-momentum approximations as for the RPA- 
meson propagators we find that the vertex is given by —i{e/gpT,.„)q'^g^^ which exactly 
corresponds to the 7p vertex in the vector dominance model of KroU, Lee and Zumino 

El- 

For the effective pvrvr coupling constant we obtain 



2 

m 



i/pTTTr ^i/Trgg ^ (o)2 ' 

This result can be used to estimate the importance of quark effects. If we take the 
commonly used value of about 6 for gp^^^^ in hadronic models we obtain for the constituent 
quark mass 

m ^ Ve/f ) . (41) 

As will be seen in the next section the leading-order pion decay constant is larger but 
not much larger than the l/A'c-corrected quantity. If we take ~ 100-150 MeV, we 
find that the constituent quark mass should be about 250-350 MeV, which is obviously 
in contradiction to our original assumption of very heavy quarks . Thus we would expect 
that this approximation does not describe the full model very well. This point will be 
discussed in more detail in the next section. 

Similar to the p-meson one can perform approximations to the l/A^^-corrected self- 
energies of the other mesons. The authors of Ref. P2[ have used effective meson-meson 



coupling constants generated from this approximation to the NJL quark loops in a linear 
sigma model calculation for in-medium pion properties. It is nice to see that a consistent 
approximation to the l/A^^c-corrected NJL model naturally generates an effective one loop 
approximation to the linear sigma model in the vr — a sector. 



A slightly different approximation to the NJL model has been performed in ref. ||15|| , 
where instead of a low-momentum expansion, the vertices are evaluated for on-shell in- 
termediate mesons. However, at least for processes dominated by intermediate pions this 
gives very similar results to those obtained with the low-momentum expansion. 
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5 Numerical results 



In this section we present numerical results for the self-energy of the p-meson and related 
quantities, such as the electromagnetic form factor of the pion and tt-tt phase shifts in the 
vector channel. Before we begin with the explicit calculation we have to come back to the 
regularization. As discussed in Sec. ||, all quark loops, i.e. the RPA polarization diagrams, 
the quark triangles and the quark box diagrams should be regularized in the same way 
in order to preserve chiral symmetry. We use a Pauli-Villars-regularization with two 
regulator masses ■,yr?7?~+~A| and -^m^ + 2A^. As before, m denotes the constituent quark 
mass and A^ is the cutoff parameter. The regularization of the meson loop (integration 
over d'^p in eq. (^TJ)) is not constrained by chiral symmetry and independent from the 
quark loop regularization. For practical reasons we choose a three-dimensional cutoff Am 
in momentum space. In order to get a well-defined result we work in the rest frame of the 
1 / A^c-iniproved meson. The same regularization scheme was already used in Ref. [17 



With the additional meson cutoff Am there are five parameters: the current quark 
mass mo, the two coupling constants gs and g^, the quark-loop cutoff Ag and the meson- 
loop cutoff Am- For a given value of Am, the current quark mass mo, the scalar coupling 
constant gs and the cutoff A^ can be fixed by fitting the pion mass m^, the pion decay 
constant /^r and the quark condensate (ipip) to their empirical values. Then we can try 
to determine the two remaining parameters, i.e. the vector coupling constant g^ and 
the meson cutoff Am, by fitting the pion electromagnetic form factor in the time-like 
region, which is related, via vector meson dominance, to the p-meson propagator (see 
below). Roughly speaking this amounts to fitting the p-meson mass and its width. The 
p- meson mass is most sensitive to the coupling constant g^. Since in the present article 
we neglect possible p and oi meson intermediate states, the value of g^ does not influence 
the observables in the pion sector, mT^^fn and {^ip)- Hence for given Am we can try 
to fit these observables together with the p-meson mass. Finally, we determine Am by 
comparing the resulting pion electromagnetic form factor with the experimental data. 

However, this procedure contains a problem: For mo 7^ the corrected pion mass m^ 
is always larger than the RPA pion mass rri§\ Hence, if we perform a fit for the l/N^.- 
improved pion mass, the masses of the RPA pions, which enter into the intermediate 
states, e.g. in diagram Fig. 0(a), are too small, i.e. the value of the threshold energy 
for the decay of the p-meson into two pions is too small. Therefore, since we are mainly 
interested in the p-meson in this article, we adjust rather than to the physical 
pion mass. Here we make use of the fact that the pion mass is very sensitive to the current 
quark mass mo whereas the two other observables we want to fit, f.^ and (ipip), depend 
only weakly on mo. 

Five parameter sets (corresponding to five different meson cutoff values Am) are listed 
in Table |l], together with the constituent quark mass m, the values of m^r, fn and {ipip) and 
the corresponding leading-order quantities. As outlined above, Ag, ttiq and gs have been 
obtained by fitting ■m^\ f^^ and (ipip) to the empirical values. For the quark condensate 
this is not known very precisely, but the absolute value is probably below 2(260 MeV)^. 
(This corresponds roughly to the upper limit extracted in ref. ||2^ from sum rules at a 
renormalization scale of 1 GeV. Recent lattice results are (ijjijj) = -2((231±4±8±6) MeV)^ 
.) It turns out that we can only stay below this limit, if the meson cutoff is not too 
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Am I MeV 


0. 


300. 


500. 


600. 


700. 


Ag / MeV 


800. 


800. 


800. 


820. 


852. 


mo / MeV 


6.13 


6.40 


6.77 


6.70 


6.54 


A 9 

9sAl 


2.90 


3.07 


3.49 


3.70 


4.16 


9v / 9s 






i.O 


i.6 


O A 

Z.4 


m 1 MeV 


ZbO. 


on /I 
304. 


396. 


44d. 


ooO. 


mV 1 MeV 


140.0 


140.0 


140.0 


140.0 


140.0 


/ MeV 


140.0 


143.8 


149.6 


153.2 


158.1 


1 MeV 


93.6 


100.6 


111.1 


117.0 


126.0 


A / MeV 


93.6 


93.1 


93.0 


93.1 


93.4 


1 MeV^ 


-2(241.1)3 


-2(249.3)3 


-2(261.2)3 


-2(271.3)3 


-2(287.2)3 


(V-^) / MeV^ 


-2(241.1)3 


-2(241.7)3 


-2(244.1)3 


-2(249.5)3 


-2(261.4)3 


(r2) 1/2 / fin 






0.59 


0.61 


0.66 




1.001 


1.007 


1.018 


1.023 


1.072 



Table 1: The model parameters (Am, Aq, mo, Qs and g^) and the resulting values of m.^, 
Jtt and (ipip) (together with the corresponding leading-order quantities), the constituent 
quark mass m and the charge radius of the pion (r^)i/^. The ratio —mo{'ip'ip)/ml.f^, is 
also given. 

large (A^^ < 700 MeV). This is related to the fact that the constituent quark mass and 
hence the absolute value of the leading-order condensate (t/^?/^) increases with Am if we 
keep fn constant. 

Table |^ also displays the ratio —mQ{'il)'ilj)/m'^fl, which would be equal to 1 if the 
GOR relation was exactly fulfilled. Note that for all parameter sets given in the table 
the deviations are less than 10% (for Am < 600 MeV even less than 3%), indicating that 
higher-order corrections in 1/Nc are small. 

We now turn to the p-meson channel. According to Eq. (|]), the polarization function 
of the p-meson is a sum of the RPA polarization loop and the four l/Ai'c-correction terms, 
as shown in Fig. |^: 

nf'"*(g) = n^'^'»^(g) + 511^^^ ^"'^''\q) . (42) 

k=a,b,c,d 

Because of vector current conservation, the polarization function has to be transverse, i.e. 

q,Ii^''^''\q) = qufl^'i^l) = 0" (43) 

For our scheme it can be shown with the help of Ward identities that these relations 
hold, if we assume that the regularization preserves this property. This is the case for 
the Pauli-Villars regularization scheme, which was employed to regularize the RPA part 
Up. Together with Lorentz covariance this leads to Eq. (|^) for the tensor structure of Up. 
On the other hand, since we use a three-dimensional sharp cutoff for the regularization 
of the meson loops, the l/A'^c-correction terms are in general not transverse. However, as 
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mentioned above, we work in the rest frame of the p-meson, i.e. q = 0. In this particular 
case Eq. (^3]) is not affected by the cutoff and the entire function lip can be written in 
the form of Eq. (^: 

fl>;'''^\q) = fl,{q)T^''6., = (U,{q) + ^ OTf (g)) T^^ , (44) 

k=a,b,c,d 

i.e. instead of evaluating all tensor components separately we only need to calculate the 
scalar functions Up = —1/3 g^ull'^'^ and ^IIp'^'' = —l/Sg^uSU^p^^'^. 

A second consequence of vector current conservation is, that the polarization function 
should vanish for = 0. For the l/A^^c-correction terms this is violated by the sharp 
cutoff. We cure this problem by performing a subtraction: 

^ SU(^\q) Yl ('^n(^^(g) - (0)) . (45) 

k=a,b,c,d k=a,b,c,d 

Note, however, that already at the RPA level a subtraction is required, although the RPA 
part is regularized by Pauli-Villars. This is due to a rather general problem which is 
discussed in detail in App. 0. 

As mentioned above, we want to fix the remaining two model parameters, i.e. and 
Am, by fitting the pion electromagnetic form factor, FTj{q), in the time-like region, which 
is dominated by the p-meson. The diagrams we include in our calculations are shown in 
Fig. ^. The two diagrams in the upper part correspond to the standard NJL description 
of the form factor if the full p-meson propagator (curly line) is replaced by the RPA 



one. Hence, the first improvement is the use of the l/A^^c-corrected p-meson propagator 
in our model. Since, in the standard scheme, the photon couples to the p-meson via 
a quark-antiquark polarization loop, in our scheme we should also take into account the 
1/iVc-corrections to the polarization diagram for consistency. This leads to the diagrams in 
the lower part of Fig. |. On the other hand the external pions are taken to be RPA pions 
(i.e. mass m^^ and pion-quark-quark coupling constant glrqq)- This is more consistent 
with the fact that the p-meson is also dressed by RPA pions and, as discussed above, we 
have fitted m^n^ to the experimental value. 

Evaluating the diagrams of Fig. ^ we obtain for the form factor 

\FA<1)\' = \ I <?/(p,p') (l - t^M)DM)) f , (46) 

where p and p' are the four- momenta of the external pions and f{p,p') is a scalar function 
appearing in the vrvrp- vertex function (see App. |C|). For the form factor it has to be 
evaluated for on-shell pions, i.e. p"^ = p'"^ = and p ■ p' = — m^. From its explicit 
form, which is given in App. 0, one can show that g^qq f{p,p') = for = 0. Since 
the p-meson self energy Hp vanishes at this point, we find -FV(O) = 1, as it should be. 

The numerical results for l-FjrP as a function of the center-of-mass energy squared are 
displayed in the left panel of Fig. ^, together with the experimental data The various 
curves correspond to different values of the meson cutoff Am- For the other parameters 
the values listed in Table |l] are taken. The vector coupling constant g^ is chosen such 
that the maximum of the form factor is at the correct position. Roughly speaking this 
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Figure 5: Contributions to the pion electromagnetic form factor. Note that the propagator 
denoted by the curly line corresponds to the 1 / N^- corrected rho-meson, while the double 
lines indicate RPA pions and sigmas. 



corresponds to fitting the mass of tlie (dressed) p-meson. We only sliow tfie results for 
Am > 500 MeV. For lower meson cutoffs the unphysical quark-antiquark decay threshold 
of the p-meson is below the maximum of the form factor and makes a comparison not 
very meaningful. 

For Km = 500 MeV (dashed line) the quark-antiquark threshold is at s = 0.63 GeV^, 
causing a cusp in the form factor slightly above the maximum. Because of the sub- 
threshold attraction in the p-meson channel the form factor drops very steeply below the 
cusp, leading to a poor description of the data above the maximum. In addition, also 
at the rising edge of the peak, where threshold effects are less important, we see that 
the width of the form factor is underestimated by the calculation, i.e. Km is too small. 
On the other hand, a meson cutoff of 700 MeV (dotted line) is already too large. In 
particular, the height of the maximum is underestimated. With Km = 600 MeV (solid 
line) we obtain the best description of the data. Since we assumed exact isospin symmetry 
in our model we can, of course, not reproduce the detailed structure of the form factor 
around 0.61 GeV^, which is due to p-cj-mixing. The high-energy part above the peak is 
somewhat underestimated, mainly due to the gg-threshold at s = 0.80 GeV^. Probably 
the fit can be somewhat improved if we take a slightly larger meson cutoff, but we are 
not interested in fine-tuning here. In addition, it is to be expected that the inclusion of 
p- and ai intermediate states will improve the high-energy behavior. 

A closely related quantity is the charge radius of the pion, which is defined as 



6 



(47) 



Results for Km = 500, 600 and 700 MeV are listed in Table |I|. All of them are close to 
the experimental value, {rlY'^ = (0.663 ± 0.006) fm p^. For Am = 700 MeV we find 



perfect agreement. On the other hand, a simple pole ansatz for the form factor leads to 
(r2)V2 = 

at q 

which are not included in our model might start to play a role. Therefore the pion charge 



p — 0.63 fm p8[, which is of the same quality as our results. Obviously, 
^ — we are not very sensitive to the details of the p-meson peak, while other effects 
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0.2 0.4 0.6 0.8 400 600 800 

s [GeV^] V s [MeV] 

Figure 6: The pion electromagnetic form factor (left panel) and the ixix-phase shifts in the 
vector-isovector channel (right panel) for different meson cutoff parameters, Am = 500 
MeV (dashed), Am = 600 MeV (solid) and Am = 700 MeV (dotted). The data points are 
taken from refs. and MW, respectively. 



radius is certainly not very well suited for fixing our model parameters. 

One can also look at the vrvr-phase shifts in the vector-isovector channel. We include 
the diagrams shown in Fig. i.e. the s-channel p-meson exchange and the direct tttt- 
scattering via a quark box diagram. The latter has to be projected onto spin and isospin 1, 
which is a standard procedure. (For example, the analogous projection onto spin and 
isospin can be found in Refs. ^.) Our results, together with the empirical data 
3Q| , are displayed in the right panel of Fig. |^. Since the main contribution comes from the 
s-channel p-meson exchange, they more or less confirm our findings for the form factor: 
below the p-meson peak the best fit of the data is obtained with Am = 600 MeV while, 
at higher energies, where gg-threshold effects start to play a role, the data are better 
described with Am = 700 MeV. 

\ ' \ / 

TT ^ / TT TT \ / TT 



TT/ ^. TT TT/ \7r 



Figure 7: Diagrams contributing to the -k-k -scattering amplitude: Quark box diagram (left) 
and s-channel p-meson exchange (right). 



Finally, we wish to compare the results of the full l/Ai'c-corrected NJL model with the 
static limit, i.e. the approximation introduced in Sec. ^. As pointed out, this approxima- 
tion corresponds to a purely hadronic description where all quark effects are suppressed. 
For instance, taking the static limit of the diagrams shown in Fig. ^, we exactly recover 
the diagrams which contribute to the pion electromagnetic form factor in the hadronic 
model of ref . . On the other hand, we already estimated that quark effects should not 
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Figure 8: The imaginary part of the 1/Nc- corrected p-meson propagator. The solid line 
represents the full calculation while the dashed line denotes the static limit. 

be negligible, i.e. we expect the static limit not to be a good approximation to the exact 
NJL calculations. In fact, if we take the parameters of our best fit to the pion electro- 
magnetic form factor {Am = 600 MeV) and insert the corresponding values for m and 
fi^^ (see Table |T]) into Eq. ( PD| ) we find gpT^Tv = 9.3, which is considerably larger than the 
value of ~ 6, typically needed in hadronic models to describe the data. Fig. || displays the 
imaginary part of the p-meson propagator for Am = 600 MeV. The solid line indicates the 
result of an exact treatment of our model. Since the parameters have been fitted to the 
pion electromagnetic form factor, the maximum is close to the empirical p-meson mass of 
770 MeV. The cusp at i/i = 892 MeV is again a gg-threshold effect. The corresponding 
result in the static limit is displayed by the dashed line. Of course there is no gg-threshold 
in this approximation. Obviously the peak is much broader and shifted to higher energies. 

For a better understanding of this behavior we also compare the real- and imaginary 
parts of the corresponding p-meson self-energies (Fig. P). Obviously, the differences seen 
in Fig. P are mainly due to the real part which is much more attractive in the exact 
NJL calculation than in the static limit, while the imaginary parts are not very different 
up to y/s ~ 800 MeV. Note that, in this region, the imaginary part exclusively results 
from the p — >• tttt process, i.e. quarks do not show up as unphysical decay channels. The 
differences to the static limit, both in the real- and in the imaginary part, come about by 
the fact, that the quarks cause a non-trivial momentum dependence of the intermediate 
meson propagators and of the effective meson-meson vertices. For instance, the quark 
triangles do not exactly behave as point-like three-meson vertices. This might be viewed 
as "physical quark effects" in contrast to the unphysical gg-decays at higher energies. 
However, the momentum dependence of any quark loop below the gg-threshold is related 
to the imaginary part above the threshold via dispersion relations. In that sense any 
difference to the static limit could be interpreted as an unphysical effect and one might ask 
whether "physical quark effects" exist at all. Of course, this fundamental question is far 
beyond the scope of our paper and can probably not be answered without understanding 
the mechanism of quark confinement itself. On the other hand, at sufficiently low energies 
(< 400 MeV) the exact self energy of the p meson does indeed almost coincide with the 
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Figure 9: The real part (left) and the imaginary part (right) of the 1/Nc- corrected p-meson 
self- energy, multiplied by —2g^. The solid line represents the full calculation, while the 
dashed line indicates the static limit. 

static limit. Here one can nicely see, how a hadronic model, which incorporates chiral 
symmetry, emerges naturally from the underlying quark structure once the many-body 
theory is carried to a sufficient degree of sophistication. 

6 Conclusions 

We have investigated meson properties within the Nambu-Jona-Lasinio model, including 
meson-loop corrections, which are generated via a systematic l/A'^c-expansion of the self 
energy in next-to-leading order. We have shown that such a scheme is consistent with 
chiral symmetry, leading to massless pions in the chiral limit. For non- vanishing current 
quark masses the pion mass is consistent with the Gell-Mann Oakes Renner relation if 
one carefully expands both sides of the relation up to next-to- leading order in l/N^. 

The relative importance of the 1/iVc-corrections is controlled by a parameter A^^, 
which cuts off the three-momenta of the meson loops. One of the main goals of the 
present article was to determine the value of Am, together with the other parameters of 
the model. To that end we have performed a fit of the quark condensate {ijjip) and the pion 
decay constant f^^, together with the pion electromagnetic form factor in the time-like 
region. The latter more or less amounts to ffiting the mass and the width of the p-meson. 
Here the meson loops are absolutely crucial in order to include the dominant p tttt- 
decay channel, while the leading-order approximation contains only unphysical gg-decay 
channels. Of course, a priori it was not clear to what extent these unphysical decay modes, 
which are an unavoidable consequence of the missing confinement mechanism in the NJL 
model, are still present in the region of the p-meson peak. 

It turns out that a reasonable fit of the above observables can be achieved with 
Am = 600 MeV. For the constituent quark mass we find m = 446 MeV. Hence the 
unphysical gg-decay channel opens at 892 MeV, about 120 MeV above the maximum of 
the p-meson peak. Our result is also interesting in the context of Ref. where we have 
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reported on the existence of instabilities in the pion sector for very large values of Aa/. 
However, with Am = 600 MeV, as determined by our parameter fit, we are far away from 
this region. In fact, we find only moderate changes in the pion and quark sector: 
and (ipip) are lowered by about 20% by the meson loop corrections, while the pion mass 
is increased by about 10%. This also indicates that our scheme converges rapidly and 
higher-order terms in the l/A^^c-expansion are small. 

We have discussed that a hadronic model can be derived from our model if one takes 
the so-called static limit. Basically, this corresponds to a low-momentum approximation 
to the quark loops. In this way one can study the importance of quark effects, which are 
present in the exact model because of lack of confinement, but absent in the static limit. 
At higher energies we find for the p-meson rather large differences between the exact 
propagator and the static limit, whereas at sufficiently low energies (< 400 MeV) the 
exact calculation and the static limit almost coincide. This implies that in the low-energy 
region the unphysical quark effects are suppressed and demonstrates that a hadronic model 
with realistic parameters indeed emerges from the underlying quark structure. It is likely 
that the quark threshold and hence the unphysical effects from continuum qq decay can 
be pushed to even higher energies once intermediate p-and Oi meson states are included. 
Given the successful vacuum description of the p meson, as presented in the present paper, 
it will be interesting to extend the calculation to finite temperature to asses medium 
modifications in the presence of a thermal heat bath. Because of spontaneous chiral 
symmetry breaking and its restoration at high temperature this implies a simultaneous 
treatment of the ai meson. Work in this direction is in progress. 
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A Definition of elementary integrals 

It is possible to reduce the expressions for the quark loops to some elementary inte- 
grals ||3^, see App. |B| and |C[ In this section we give the definitions of these integrals. 

J _ [ d^k 1 

^^^^ ^ j (27r)4 (F - m2 + ie){{k + pf - + le) ' ^^^^ 
^^^^ ^ j (27r)4 - m2 + iey{{k + pf - + is) ' 

M{pi,p2) = j (27r)4 (F -m? + ie){kl - w? + ze){kj - w? + le) ' ^^^^ 
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J, [ d^k 1 

(Pi,P2,P3) - j (2vr)4 (p - m2 + ie)(A;2 - m2 + ie){kl - + - m2 + ie) ' ^^^^ 

(2^ (P - m2 + (fc? _ + _ + ^ (53) 

with ki = k + pi. The function Mi{pi,p2) can be expressed in terms of the other integrals: 
, . P1-P2 I{pi) -pI I{P2) + {pl-Pi-P2) I{pi-P2)+pl {pl-PrP2) M{pi,p2) 

Ml[Pl,P2) — r-77 2~^T , 

2 {{prP2y - Ptpi) 

(54) 

All integrals in Eqs. (^8|) to (|53D , are understood to be regularized. In our model we use 



Pauli-Villars regularization with two regulators, i.e. we replace 



i=o 



with 



/i| = + j Ag ; Co = 1, Ci = -2, C2 = 1 . (56) 

One then gets the following relatively simple analytic expressions for the integrals Ji, I{p) 
and K{p): 



..,^^E.(..K.,...M-..,..M-..O...M.S.)) ... 

j 

A^^ = 0) = I^E^^1^(/^|) (59) 



^^^^ " 16^ ^ 2p2(a;^.^^_a;.^) ( " ^^(^^^^ ~ ^'^(-^Ji) + ^^(^^2) + ln(-x,-2)) , (60) 



with 



1 + Ml 



a;,i,2 = -±-Wl-^. (61) 



An analytic expression for the three-point function (Eq. ^T]) can be found in Refs. 
and In certain kinematical regions the four-point function (eq. ^2]) is also known 



analytically |33, |34 



B RPA propagators 

Using the definitions given in the previous section the gap equation (Eq. (|^)) takes the 
form 

m = mo + 2igs AN^Nf m h . (62) 
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Similarly one can evaluate the quark- antiquark polarization diagrams (Eq. (^) and cal- 
culate the RPA meson propagators. The results read: 

^'^^^^ ^ ^ + 2z^7. 2iV,iv/(p2_4m2) I{p) ' ^^^^ 

^'^^^^ ^ ^ + 2zg, 2N^Nf f Hp) ' ^^^^ 

n = fRc^N 

''^^^ 1 + 2ig, iN.Nf (-2m2 /(O) + (p^ + 2m2) I{p)) ' ^ ^ 

(j)) = (QQ) 

"'^^^ 1 + 2tg, iNcNf (-2m2 /(O) + (p^ - 4m^) I{p)) ' ^ ^ 

We should comment on the p- and the ai-propagator. A straight-forward evaluation of 
the vector polarization diagrams gives 

U,{p) = -t^N.Nf (-2/1 + + 2m^) I{p)) . (67) 

Because of vector current conservation this function should vanish for p'^ = 0. This is only 
true if 

/(O) = /i , (68) 

which is not the case if we regularize I{p) and Ji as described in App. 0. This corresponds 
to the standard form of Pauli-Villars regularization in the NJL model 0. Alternatively 
one could perform the replacement Eq. ( |55D for the entire polarization loop. In fact, this 
is more in the original sense of Pauli-Villars regularization ||35|. Then the factor in 
Eq. (^) should be replaced by a factor inside the sum over regulators and one can 
easily show that Eq. ( p8| ) holds (see Eqs. (|57D and (|^)). However, this scheme would 
lead to even more severe problems: From the gap equation (Eq. |6^) we conclude that ili 
should be positive. On the other hand the pion decay constant in the chiral limit and in 
leading order in l/N^ is given by 

/2(o) = -2iNcNf J(0) . (69) 

which implies that im^li^^ should be negative. So irrespective of the regularization 
scheme Eq. (BSl) cannot be fulfilled if we want to get reasonable results for m and at 



the same time. Therefore we choose the standard form of Pauli-Villars regularization in 
the NJL-model [|^ and replace the term Ji in Eq. (^) by hand by rn? /(O). This leads to 
the p-meson propagator as given in Eq. (|65|) . For consistency the a\ has been treated in 
the analogous way. 



C Explicit expressions for the meson-meson vertices 

In this section we list the explicit formulae for the meson-meson vertices. We restrict 
ourselves to those combinations which are needed for the calculations presented in this 
article. 
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We begin with the three-meson vertices TMi,M2,M3ilyP) (^^e Fig. H): 
-ir,,,,,(g,p) = i2mN(^I{p') + I{q) + I{p) + {4m^-^{p'^ + p^ + q^))M{p,-q)^ , 

-tTl'^^,iq,p) = i2mN5a,[l{p')+p-qM{p,-q)Y 

■-rMf ^ X \f P^q^q^ + g^pV - p-q{p^q^ + q^p^)\ 
-2r;^,,;,(g,p) = Sa,Kq,p)[g^ p2g2 _ (p.g)2 j' 

h{q,p) = imN(^I{q) + I{p)-2I{p') + {Am^ -2p-q-p'^)M{p,-q)^ , 

-ir^;S(g,p) = ^a6c(g^/(g,p) -p'^/Ip,?)) , 

= iv(-/(g)+p2M(p,-g) + 2p-gMi(g,-p)) , (70) 

with p' = —p — q and = AN^Nf. 

For the four-meson vertices we only need to consider the special cases needed for the 
diagrams (b) and (c) in Fig. 

-^r.,.,.,.(g,p,-g) = -n[ ~ + + + Am' {M{p, q) + M(p, -q)) 

+2(m^ (4m^ — p^ — q') — ) — g, p — g) | 

-2r,,,,,,,(g,p,-p) = -Ar|/(p + g) + /(0)+4m2(K(p) + ir(g) + 2M(p,-g)) 

+2p-qM{p, -q) - q'K{q) - 

2 2 

+m2(16m2 - V - V + ^)L{P, -q, 0)| 
-2r:j|'^^,^^(g,p,-g) = (5,bAr|/(p + g) + /(p-g)+p2(4m2-g2)^(p^_g^p_g)| 

-tTl'^^^^^{q,p, -p) = 5atN{ - lip + g) - /(O) - (4m^ - q'){K{q) - p'L{p, -g, 0)) 

+p'K{p) - 2p-q M{p, -g)| 

-zr^^i_^(g,p,-p) = -Nnabcd{l{p + q) + nO)-p'Kip) 

-q'K{q) + 2p-qM{p, -q) + p'^q^L^p, -q, 0) | 

-^r^;',_^_,(g,p, -g) = -25„biv|/(p + g) + /(p - g) + 2/(g) - p-g(M(p, -g) - M(p, g)) 

+ (4m^ - 2p^)(M(p, g) + M(p, -g)) 
W(8m2 _ g^2 ^ 4^2 ^ P_z(P:3L^L{p, -q,p- g)| 

-^r^^,,,^^(g,p,-g) = -26abN[-Iip + q)-I{0) + {p'-4m')K{p) 

+ (g^ + 2m')K{q) + {Am' - 2p-q)M{p, -g) 
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W(8m' - 2/ + V - ^)L{P, -q, 0)1 

-iTf;'^^^{q,p, -q) = 2NKatcd{ - I{p + q) - I{p - q) - 21{q) 

+2p'^{M{p, q) + M(p, -q)) + p-q{M{p, -q) - M{p, q)) 
+ (2m^p^ -p'^ + {p-qY)L{p, -q,p - g)| 

-tT';^'^%^^{q,p,-q) = 2NKabcd{lip + q) + IiO)-p'Kip)~{q' + 2m')Kiq) 

+2p-qM{p, -q) + p\2m^ + q^)L{p, -q,0)^ , (71) 

with rp,M,M,p(g,p, -g) = ^;.i.r;^^M,M,p(?>p> -?)>rp,M,p,M(g,p, -p) = 9f^>yT''p'^M,p,M(i^P^~p'> 
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